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Abstract 

A solution to the Schrodinger equation for the nonrelativistic Green function which 
is used for describing the heavy quark-antiquark pair production near the threshold in 
e^e~ annihilation is presented. A quick comparison with existing results is given. A 
choice of the effective mass scale for the nonrelativistic system with Coulomb interac- 
tion is discussed. 
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The hadron production near heavy quark threshold will be thoroughly studied experi- 
mentally at future accelerators, e.g. [Q. The dynamics of a slow moving pair of the heavy 
quark and antiquark near the production threshold is nonrelativistic to high accuracy that 
justifies the use of the nonrelativistic quantum mechanics as a proper theoretical framework 
for describing such a system [0]. Being much simpler than the comprehensive relativistic 
treatment this approach allows one to take into account exactly such essential features of 
the dynamics as Coulomb interaction ^ . The spectrum of hadronic states produced near 
the quark-antiquark threshold is contained in the Green function G{E) = [H — E)^^ of the 
effective nonrelativistic Hamiltonian H. In the position space it reads 

G{E-vy) = {v\{H-E)-y) (1) 

and satisfies the Schrodinger equation 

{H -E)G{E-T,r') = 5(r-r'). (2) 

The Hamiltonian H is represented in the following general form H = Hq + AH where the 
first term Hq is the Coulomb Hamiltonian 

H, = P--^ (3) 

m r 

with p'^ = — A^. The color factor for the fermion representation of the gauge group is 
Cp = {N^ — l)/2Nc, as is a QCD coupling constant, m is a mass of the heavy quark. The 
normalization point for as will be fixed later. The second term AH accounts for relativistic 
and perturbative strong interaction contributions which are assumed to be small, i.e. they 
are treated as corrections to the Coulomb spectrum. The term AH has an explicit form 

AH = AkV + A.potV + AnaV + AbfV. (4) 

The quantity AkV in eq. (HI) is a relativistic kinetic energy correction 



the second term ApotV is the strong interaction perturbative corrections to the Coulomb 
potential 0, ^, the third term AbfV is a Breit-Fermi potential with addition of the color 
factor Cp 0, §]. AnaV = —CACFal/{2mr^) is the non-Abelian potential of quark- ant iquark 
interaction |Q, Ca = is a. color factor for the gluons. The corrections to the Coulomb 
Green function at the origin due to terms A^V, V^a and Vbf have first been presented in 
1^, [Tl[|. Numerous applications of these results and further references can be found in recent 



literature, e.g. [|T2|, |T^, |15|, The treatment of the problem in the present paper is 
slightly different. For the s-wave production of a quark-antiquark pair in the triplet spin 
state {L = 0, S = 1) the Breit-Fermi potential reads |0, |^ 



The correction AH can be be rewritten as 



AH = ApotV - -^Hl 
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where ip^ = 9^ + 1/r. Note that a similar form of the representation of the correction in 
terms of powers of the leading order operator is usual for anharmonic oscillator problems 



\% . The part of the potential proportional to the 5-function 5{r) is a separable potential. 
The equation for the full GF with such a potential can be exactly solved for the quantity we 
need. Noticing this fact one rewrites the Hamiltonian in the form 

H = Ho + AH = Hir + asV (7) 

with 

Here ifj^ is an irreducible Hamiltonian. The new representation for the GF reads 

G{E) = {H - E)-^ = {Hir + a^V - E)-K (9) 



Introducing the GF Gir{E) of the irreducible Hamiltonian 

G,r{E) = {H,,-E)-^ (10) 

one obtains the following equation for the full Green function G{E) (written in the operator 
form) 

G{E) = Gir{E) - asG„{E)VG{E) . (11) 
Eq. ( [Til) is exactly solved for the position representation component G{E] 0, 0) with the result 

1 + asVoGir{E;0,0} 

Eq. (|T2|) accounts for the 6{r) part of the correction to the irreducible Hamiltonian exactly. 
It is analogous to the usual representation of the vacuum polarization function through the 
one particle irreducible block. In this case the irreducible object Gir{E; 0, 0) is characterized 
by the absence of 6{r) interaction in it. The function Gir{E; 0, 0) can be found perturbatively 
using the Coulomb solution as a leading order approximation. 

The Green function G{E) emerges as a nonrelativistic limit of the relativistic scatter- 
ing amplitude near the production threshold. The nonrelativistic Hamiltonian can be con- 
structed from the QCD Lagrangian. In the leading order of nonrelativistic expansion there 
is an energy independent factor (matching coefficient) G{as,m) that allows one to map the 
quantum mechanical quantities onto the relativistic cross section near the threshold. In 
higher orders of nonrelativistic expansion further terms of the expansion of the current itself 
and new vertices of the effective Lagrangian are generated that should be accounted for when 
the cross section in QCD is calculated. We do not discuss these terms here because their 
contributions start at the higher orders in expansion parameters. Therefore the correspond- 
ing expressions should be taken only in the leading order of hard loop expansion. The cross 
section of the heavy quark-antiquark pair production near the threshold in e~^e~ annihilation 
contains a part with the nontrivial loop expansion. It has the form 

R^'is) ^ G{as,m)lmG{E;0,0) , s = {2m + E)\ (13) 



Here C(as, m) is the matching (hard or high energy) coefficient, ois is the couphng constant. 
Generally, one can use the different normalization points (or different subtraction proce- 
dures) for Us and for the corresponding coupling constant which enters the expression 
for the nonrelativistic Green function. The vacuum polarization function near the thresh- 
old C {a siTn)G{E\ requires subtraction. This feature is familiar from the PT analysis 
of the vacuum polarization function in the full theory. To obtain a finite quantity one 
can differentiate C(as, m)G(i?; 0, 0) with regard to E (constructing the D-function) or take 
the discontinuity across the physical cut (constructing the imaginary part or cross section 
R^^{s)). Both D-function and R^^'{s) are finite. We write 

and obtain the following expression for the D-function 

C(.„ ,n)±G(E, 0, 0) ^ Ag..(B; 0, 0) . (14) 

The imaginary part of C(as, m)G{E; 0, 0) reads 

n^- M nrj? n n\ C(a„m)Im Gir{E;0,0) 

C(a,. ™)I,„ G(Bi 0. 0) = ,i+„,v',ReG„(J5iO.O))= + (a,V„ImG4B;0.0))^ ' ^'"^ 

The quantities in eqs. (0) and (|15D are finite. The explicit expression for the coefficient 
C(a<j,m) has been found at the a"^ order within dimensional regularization ||18|. It contains 
a singularity of the form 

Gsias, m) = 1 - Gpal + ^) 1 (16) 

where only the singular part Gsi^s, ^) of the coefficient G{as, rn) is written. This singularity 
cancels in eqs. (0) and ([T5|). We consider this cancellation (the renormalization procedure) 



for the case of D-function only. To fulffil the renormalization procedure we use the dimen- 
sionally regularized Green function in the expression for the correction proportional to Vq. It 
suffices to substitute the pure Coulomb Green function for this purpose. The dimensionally 



regularized Coulomb Green function at the origin takes the explicit form |]T5 



0. 0) ^ - + ^ ^- + .„ ^ - 2* ^1 - j J I (17) 

where = —mE, ip{z) = r'{z) /r{z) is digamma function and T{z) is Euler's F-function. 
One finds 

(1 + asVoGir{E; 0, 0))^ ^ (1 + a, 0, 0))^ = Z'^l + asVoGc{E; 0, 0))^ (18) 

where 

Z-' = l + C^m'n- = 1 - Cal + ^) - . (19) 
An e \ 6 2 / e 

The finite (renormalized) Coulomb Green function at the origin has the form 

,,,.;0.0)^£{-.. ^(.4-2,(1-^))}. ,20) 

The renormalization constant Z cancels the divergence C5(a<j,m) of the high energy coeffi- 
cient C{as, in a proper order of loop expansion. One uses as = as + O(a^) as a formal 
PT relation to achieve the cancellation. The finite coefficient of order (or a^) depends 
on particular ways of subtraction in C{as, ^) and Gir{E). If the subtraction procedures for 
C{as, ITT') and Gir{E) are not properly coordinated during the calculation the finite coefficient 
is fixed by matching [|Ty|. If the calculation for both C(a^5,m) and Gir{E) has been done 
within one and the same subtraction scheme this matching is automatic, e.g. [^. Therefore 
after the standard renormalization, eq. (|T2|) gives the representation of G{E; 0, 0) as a Dyson 
sum of irreducible terms. The spectrum of the full system is determined by the equation 

Gir{E-0,0)-^ + asVo = (21) 

where Gir{E; 0, 0) is constructed perturbatively. Eq. ( 121|) can be solved exactly or perturba- 
tively. The isolated roots of eq. (|2T|) give the discrete spectrum of the system. There is also 
a continuous spectrum given by the discontinuity of Gir{E;0,0) across the cut at positive 
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values of E. For the continuous spectrum one can use the Coulomb GF in the denominator 
of eq. (^). 

The problem of calculating the near-threshold cross section reduces to the construction 
of the spectrum of the irreducible Hamiltonian Hi^ 



Hir = Hq + iSVpot — -f^o ~ 



r 



4m 4m 

The spectrum is found within PT. The leading order spectrum is given by the renormalized 
pure Coulomb solution eq. (^). The term ApotV represents the first and second order 
perturbative QCD corrections to the Coulomb potential which were studied. The correction 
due to the first iteration of the ApotV term has been found in ref. |^ where the simple 
and efficient framework for computing the iterations of higher orders was formulated. The 
explicit formulas for these corrections can be found in |]T6| . For the kinetic Hq term one finds 



Pk{E) = Y^5iE'-^~E)piE') (23) 



where p{E) is a density of Coulomb states with the energy E 



p{E) = -lmGc{E;0,0). (24) 

TT 



The sum is over the whole spectrum. For the discrete levels it gives 



Pk{E)Usc = E'^(^n - ^ - E)\MO)\^ (25) 

4m 



n 



with ijjn being a bound state with the energy level E^- The position of the pole is now at 
EP"''^ = En — E^/Am. For the continuous spectrum the correction reads 



E'^ _ 1 



Pkmcont = / dE'6iE' - — - E)p{E') = ^Tn^piE) (26) 

Jo 4m 1 — E / 2m 

with 

E- E + = . (27) 

4m- 
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Eq. ( p7D can be solved perturbatively for small E. One obtains 

E^ 



E = E + — + O(E^) (28) 



which is the substitution used in refs. [|^, 0. Note that there is no correction to wave 
functions because has no non-diagonal matrix elements between the Coulomb states. The 
last term in eq. (|22|) gives no correction to the spectrum. The term with the anticommutator 
[HQ,r~^]^ in eq. (|22|) gives the correction to the spectrum which was presented earlier as an 
energy dependent shift of the coupling constant. This correction can be written in the form 

Gc{E) + AaGc{E) = Gc (e; a, (^l + (29) 

which is valid at small E and was presented in |T^, Note, however, that while the 
^-function part of the correction to the Hamiltonian can be considered unique because it 
represents a reducible vertex, the irreducible potential Hir can be chosen in different forms. 
Indeed, one can rewrite the sum of kinetic and anticommutator terms of Hir in the form 



1 ^-^o 3Cpfto 
-—Hi 



Am ^ Am 



Ho, - 
r 



5 „^ 3 



——H 

Am ^ Am'^ 



Ho,p' , (30) 



which is a reshuffling of contributions within the irreducible correction. While the kinetic 
term Hq has only its coefficient changed as a result of such a reshuffiing, the new anticom- 
mutator term leads to a modification of the leading order GF of the following form 

1 3 / 2 1 1 2\ 1 ^ 3p2 1 
+ —-J ip^— + — -p^ + — -E- ^ 



Ho-E Am^ \ Ho-E Ho-E"^ ) Ho-E 2m^ Ho - E 
3 / 2 1 1 2\ (P^ f, 3E\ Cpas ' 



The first term of this equation does not affect the structure of the spectrum. The last 
term in eq. (^) can be interpreted as a correction to the mass. To a certain degree this is 
equivalent to the previous case (the correction to the coupling in eq. (^)) because the 
genuine parameter of the Coulomb problem is the Bohr radius or momentum, Pb = oigm (or 



CFOis'm/2). One can see that both changes as — (1 + 3E/2m) and m — > m/ (1 — 3E/2m) 
lead to the same result for ps within the accuracy of the approximation, i.e. up to higher 
order terms in E 

PB-P^(l + g)=YZ^ + 0(i?^). (32) 

One can check that the total correction to the discrete energy levels, for instance, is the 
same. The old solution eq. ( pOf ) gives 



mAE„ = ~El + 3El = ^El (33) 



while the new decomposition eq. (|3TD results in 



mAE^ = lEl + lEl=^-lEl. (34) 

This coincidence is valid only parametrically in the region of applicability of perturbation 
theory in the expansion parameters E/m and/or as- It can result in numerical difference 
when extrapolated to larger energies in the continuous spectrum. The different forms of 
the decomposition may lead to different numerical predictions within numerical evaluation 

The expression given in eq. (|T2|) has a standard structure of the polarization function. 
The solution to eq. (^) (one or the first of the poles of the generalized "propagator" of the 
full system with the Hamiltonian H) is an important dimensional parameter for the system. 
For the observables which are saturated with the contributions of the discrete spectrum the 
quantity Ef can serve as a natural mass parameter. In this case the position of a pole (the 



numerical value of some solution Ef to eq. (|2l|) ) can be chosen as a scale for the system 



instead of the heavy quark pole mass. However, for the observables which are saturated 
with the contributions of the continuous spectrum or have a considerable admixture of such 
contributions, the natural mass scale is not necessarily related to Ej and is determined by 
other properties of the full interaction. 
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To conclude, we have presented a solution to the Schrodinger equation for the nonrel- 
ativistic Green function. The solution has the form of a resumed geometric series (Dyson 
resummation) of irreducible blocks which is usual for quantum field theory. The property 
of irreducibility is defined with respect to the 5-function part of the interaction potential. 
Different decompositions of the irreducible Hamiltonian Hir for the treatment within PT 
are considered. A new form of the first order correction to the spectrum of irreducible 
Hamiltonian Hir is given. 
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